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Predicting  the  Liquid  Lengths  of  Heavy  Hydrogen  Fuels 


Laura  L.  Hoogterp 


Abstract 

The  purpose  of  this  paper  is  to  outline  the  procedure  used  in  determining  the 
liquid  tmgths  in  diesel  fuels.  Using  models  formulated  by  previous  researchers  as  well 
as  the  tliermodynamic  properties  for  three  fuel  surrogates  the  liquid  length  can  be 
determined  for  diesel  fuel,  JP8  as  well  as  provide  a  model  for  fuels  between  these  ranges. 
This  information  is  to  later  be  used  to  develop  a  computer  program  to  perform  these 
calc:ulations  on  engines  running  in  a  test  cell. 


Introduction 

In  order  to  predict  the  liquid  length  in  diesel  fuels  three  surrogate  fuels  are  chose 
to  model.  Dodecane,  C12H26,  is  used  to  model  JP8,  Hexadecane,  C16H34,  models 
diesel  fuel,  and  Tetradecane,  C14H30,  is  used  to  model  fuel  between  these  ranges. 

Many  people  have  previously  done  research  into  predicting  the  liquid  lengths  of 
diesel  fuels.  This  information  is  then  used  to  help  predict  liquid  lengths  for  the  three 
fuels.  The  models  used  are  taken  form  the  SAE  technical  paper  by  Dennis  L.  Siebers 
while  the  thermodynamic  information  on  the  fuels  is  found  in  books  published  by  the 
American  Petroleum  Institute  and  a  thermodynamics  book. 

These  sources  are  used  to  find  the  liquid  lengths  which  are  then  checked  against 
existing  data  for  hexadecane.  Once  the  process  is  determined  to  be  valid,  the  same 
procedure  is  followed  for  both  tetradecane  and  dodecane.  The  method  used  is  then  broke 
down  so  that  it  can  be  used  in  a  computer  program  that  will  predict  the  liquid  lengths. 

The  program  will  be  used  on  engines  that  are  running  in  a  test  cell.  The  rest  of  this  paper 
will  outline  the  procedure  that  is  used  to  predict  the  liquid  lengths  in  these  fuels. 

Modelling  the  Liquid  Length 

In  order  to  determine  the  liquid  length  of  fuels  equations  that  have  already  been 
formulated  through  previous  research  are  used.  According  to  the  SAE  paper  by  Siebers 
the  following  equation  is  used  to  solve  for  the  liquid  lengths  [1]: 


-1 


whe:re 

L  is  the  liquid  length  of  the  fuel, 
b  is  a  constant  stated  to  be  0.41, 
a  is  a  constant  found  to  be  0.66. 

Pj^  is  th  e  density  of  the  liquid  fuel, 

is  the  density  of  the  ambient  gas. 


Q  is  the  area  contraction  coefficient, 

d  is  the  diameter  of  the  orifice, 

9  is  the  spray  spreading  angle,  and 

B{T^ ,F\,Tf  )  is  the  gas  mass  flow  rate  based  on  both  enthalpies  and  compressibilities  and 
sobbed  for  iteratively. 


This  value  of  B  is  solved  for  using  the  following  equation  [1]: 


Z,iT,.P,HP,-P.]-M^  h,(T,)-h,(T,,P,) 


where 

<Ts  5  )  is  the  ambient  air  compressibility  based  on  saturation  temperature  and 

partial  pressure, 

Zj  {T^ ,  Pj )  is  the  vaporized  fuel  compressibility  based  on  the  saturation  temperature  and 
pressure, 

is  thi;  saturation  pressure  of  the  fuel, 

Mj.  is  the  molecular  weight  of  the  fuel, 

[P„-P,]  is  the  partial  pressure  of  the  ambient  gas,  and 
is  the  molecular  weight  of  the  ambient  gas. 

On  the  (mthalpy  side  of  the  equation, 

h^iT^,P„)  is  the  ambient  gas  enthalpy  based  on  initial  ambient  temperature  and  pressure, 
K  (-^s  ’  “  ^5 )  is  ambient  gas  enthalpy  based  on  saturation  temperature  and  partial 

pressur(5, 

hj^{T^)  is  the  fuel  enthalpy  based  on  saturation  temperature,  and 
hj{Tf,P^)  is  the  fuel  enthalpy  based  on  initial  fuel  temperature  and  pressure. 

The  compressibility  of  the  ambient  gas  is  assumed  to  be  one.  This  assumption 
can  be  made  when  the  pressure  of  the  gas  is  low  when  compared  to  the  critical  pressure 
and/or  the  temperature  is  large  relative  to  the  critical  temperature  [3].  One  or  both  of 
these  conditions  are  met,  which  validates  the  assumption. 

The  compressibility  of  the  fuel  when  it  is  a  saturated  vapor  is  found  by 
interpolating  numbers  off  a  table  found  on  the  American  petroleum  institute’s  technical 
data  book.  To  simplify  this  procedure  so  that  it  can  be  transformed  into  computer  code  a 
best  fit  curve  of  the  data  is  found.  A  plot  is  created  of  the  compressibility  factor  based 
on  the  reduced  saturation  temperature.  A  polynomial  curve  fit  is  then  performed  to 
obtain  a  relationship  between  temperature  and  compressibility  for  each  fuel.  These 
relationships  should  not  be  used  outside  the  temperature  range  of  450  K  to  1300  K  to 
ensure  the  integrity  of  the  curves. 

For  dodecane  the  relationship  used  to  determine  the  compressibility  is 
=-16.85-7;',  +36.104-r^',  -  26.425  •  ,  +  7.5406 

where  ,  is  the  reduced  saturation  temperature  of  the  fuel. 

The  tetradecane  compressibility  relationship  is 
Z^  = -1 7.924 -r,',  +36.143-7;',  -  24.71 -T,  ,  +6.6857 

and  hexadecane’s  relationship  is 

Z^  =  -16.587  -  T;',  +  34.594 -  -  24.531  - 7; ,  +  6.869  . 


The  enthalpy  of  the  ambient  gas  is  looked  up  in  tables  found  in  reference  3  and 
are  based  on  either  initial  or  saturation  temperatures.  To  make  calculations  easier  a 
relationship  is  formulated  to  determine  these  values.  Since  the  enthalpies  are  taken  over 
a  wide  range  of  temperatures  they  are  broke  up  into  smaller  sections  to  keep  the  equation 
more  accurate,  and  should  not  be  used  outside  the  stated  temperature  ranges.  The 
enthalpy  is  in  kJ/kg  and  the  temperatures  are  in  Kelvin. 

For  temperatures  ranging  for  450  K  to  600  K  the  enthalpy  is  found  using 
=:  1.033*7-13.605. 

For  temperatures  ranging  from  600  K  to  800  K 
h„  ==1.0735*7-37.788, 

temperatures  between  800  K  and  1000  K  use 
==1.1184*7-73.254, 


and  from  1000  K  to  1320  K  use  the  relationship 
==1.1682*7-123.49. 


To  determine  the  enthalpy  of  the  fuel  in  a  gaseous  state,  a  procedure  is  used  out  of 
the  American  Petroleum  Institute  data  books.  Two  double  integrations  are  performed  and 
then  the  results  plugged  into  an  equation.  Due  to  the  complexity  of  this  procedure  the 
results  are  also  plotted  with  a  best  fit  curve  for  fuel  in  a  gas  and  liquid  state.  For  fuels  at 
a  reduc(2d  pressure  less  than  0.2  and  in  a  gas  state  the  following  equations  can  be  used 
[2]: 

dodecane;  =  1562.5  ■  7^  ^  -  444.84 

tetradecane:  hj-^  =  1683.8  •  7^  ^  -  469.88 

hexadecane:  =1869.7 -7^5 -550.15 

where  ^  is  the  enthalpy  of  the  fuel  when  it  is  a  vapor. 

If  the  reduced  pressure  of  the  vapor  is  greater  than  0.2  then  the  following 
equations  listed  below  are  used.  For  the  set  of  conditions  given  this  occurs  at  550  K  in 
dodecane,  580  K  in  tetradecane,  and  620  K  in  hexadecane.  The  equations  are: 

dodecane:  h^-^,  =-\92\.9-T^^  +4996.7-7^^  -1978 
tetradecane:  = -3712.2 -7^^^  +  8085.6- 7^  ^  -3229.2 

hexadecane:  = -929.51 -7^^^  +3520.3-7^^  -1283.4 

where  h  is  given  in  kJ/kg. 


To  deteimine  the  enthalpy  of  the  fuel  in  a  liquid  state  the  following  equations  are  used. 
They  are  based  on  pressure  rather  than  temperature  since  the  fuel  entering  the  system  is  at 
a  constant  temperature.  Enthalpy  is  in  kJ/kg. 

dodecane:  =  1.277  •  ^  +  316.3 

tetradecane:  hj-i  =1.1718-P^^ +319.12 
hexadecane:  =  1.0429  •  +  320.95 

where, 

^  is  tlie  reduced  ambient  pressure  with  respect  to  the  fuel  and 
hj  I  is  tlie  enthalpy  of  the  fuel  in  a  liquid  state. 

To  calculate  the  spray  spreading  angle  the  following  equation  is  used  [1]: 


=  c- 

^  Pa^ 

0.19 

-.0043- 

h 

.2, 

[Pf] 

V 

p. 

where 

c  is  a  constant  based  on  the  orifice  diameter  and  is  0.26  for  a  diameter  of  246  fim. 
To  determine  the  fuel  density  in  its  liquid  state  the  following  equation  is  used  [2]: 


where 

is  the  saturated  liquid  density, 

P^  is  the  critical  pressure  of  the  fuel, 

Tj- ,  is  the  fuel  temperature  reduced,  and 

is  an  empirically  derived  constant  that  can  be  found  in  reference  2. 

Since  the  equation  came  from  reference  2  and  the  formulation  of  the  constants  is 
unknovm,  the  critical  pressure  must  be  inputted  in  psi.  This  produces  density  with  units 

of  and  can  then  be  converted  to  SI  units. 

The  ambient  and  saturation  pressure  are  used  in  determining  many  properties.  In 
order  to  find  the  ambient  pressure  before  fuel  enters  the  system  the  ideal  gas  law  is  used 

[3]: 


P.=PaRT„ 


where 

R  is  the  ideal  gas  constant, 
is  the  ambient  pressure,  and 

is  th(e  ambient  temperature. 

To  ilnd  the  vapor  pressure  of  the  fuel  based  on  a  its  temperature  the  correlation 
In P  =  .4  +  —  +  ClnJ;  +  Dr/  +  ^  is  used  [2], 

3  rri  J  A  *■ 

^  s  ^  S 

A,  B,  C,  D,  and  E  are  constants  based  on  the  fuel  in  question.  It  should  be  noted  that  the 
T  should  be  in  degrees  Rankine  and  the  equation  produces  a  pressure  in  psi. 

Results 

A  spreadsheet  in  excel  is  set  up  in  order  to  perform  the  calculations  necessary  to 
dete;rmine  the  liquid  length  of  the  three  fuels.  In  order  to  check  the  results  from  excel 
with  the  work  previously  done  by  Siebers  the  same  set  of  conditions  are  used.  The 
ambient  gas  densities  used  are  3.6,  7.3,  14.8,  30.2,  59.0  kg/m^. 

In  Seiber’s  research  the  ambient  gas  conditions  have  a  composition  of  89.7%  N2, 
6.5%  CO2,  and  3.8%  H2O  with  a  molecular  weight  of  28.67  [1].  The  excel  sheet  uses  air 
as  its  ambient  gas  and  the  appropriate  thermodynamic  properties.  This  difference  can 
cause  a  difference  of  4%  in  the  final  liquid  length  due  to  heat  capacities  [1]. 

The  injected  fuel  temperature  is  assumed  to  be  438  K,  and  the  orifice  diameter  is 
246  jitm.  The  liquid  length  is  then  measured  at  700  K,  850  K,  1000  K,  1150  K,  and  1300 
K.  Using  a  table  in  Siebers  report  the  area  contraction  coefficient  is  assumed  to  be 
0.8124. 

These  conditions  coupled  with  the  equations  listed  above  are  used  to  iteratively 
soh^e  for  B.  The  conditions  are  determined  at  a  variety  of  saturation  temperatures  at  the 
each  of  the  ambient  temperatures  and  densities.  The  equation  for  B  is  then  solved  to 
witliin  a  hundredth  of  each  other.  It  is  at  this  point  that  the  liquid  length  can  be 
determined.  Table  1  is  a  sample  of  the  excel  sheet  created. 


Table  1 :  Excel  sheet  created  to  find  the  liquid  length  of  hexadecane  with  ambient 


conditions  at  850  K  and  a  density  of  30.2  kg/m^ 


Knowns: 

f^air 

fkq/m^31 

Pressure 

a 

b 

Ca 

Mair 

Mfuel 

d 

tan  9/2 

Pa  [Pa] 

Pa  [Mpa] 

0.66 

0.41 

30.2 

0.8124 

28.97 

226.45 

0.00025 

0.139337 

7366944 

7.366944 

/sfuel 

Tc.air 

Tc.fuel 

Pc, air 

Pc, fuel 

Pr, ambient 

accentric 

Tf[K] 

Ta[K] 

fkq/m''31 

[K] 

[K] 

[Mpa] 

[Mpa] 

Tr.fuel 

[Mpa] 

factor 

43ii 

350 

662.939 

133 

723 

3.77 

1,40006 

0.607864 

5.2618848 

0.7174 

literal  ons 

for  “s 

Tr.fuel 

Ps 

Pr, fluid 

Pa-Ps  Za(Ts,Pa-Ps)Zf{Ts,Ps)  ha(Ta,Pa)  ia(Ts,Pa-Ps  hf(Tf,Pa) 

hf(Ts) 

B(z) 

B(h)  difference  Average  B  Liquid  length 

[Kl 

fl 

[MPal 

[1 

fMPal 

f1 

fl 

mm 

[kJ/kq] 

[kJ/kql 

[kJ/kq] 

[1 

11 

fl 

fl 

[ml 

450 

C.622 

0.00427 

0.0031 

7.363 

1 

0.97297 

877.175 

451.8 

325.74 

626.4521 

0.004664 

1.4 

-1.4099 

0.709612 

0.01706802 

500 

C'.692 

0.02256 

0.0161 

7.344 

1 

0.95738 

877.175 

503.02 

325.74 

744.8596 

0,025077 

0.9 

-0.86764 

0.458897 

0.02438106 

550 

C'.761 

0.08132 

0.0581 

7.286 

1 

0,92846 

m.m 

554.74 

325.74 

870.8941 

0,093975 

0.6 

-0.49748 

0.342716 

0.03131938 

600 

0.83 

0.22607 

0.1615 

7.141 

1 

0.85258 

snm 

607.02 

325.74 

1002.705 

0.290263 

0.4 

-0.1088 

0.344665 

0.03116484 

650 

Ci.899 

0.52573 

0.3755 

6.841 

1 

0.72489 

877.175 

659.84 

325.74 

1123.17 

0.82867 

0.3 

0.556126 

0.550607 

0.02094906 

700 

0.968 

1.08217 

0.7729 

6.285 

1 

0.48614 

877.175 

713.27 

325.74 

1251.834 

2.768636 

0.2 

2,591651 

1.472811 

0.00989018 

610 

Ci.844 

0.27092 

0.1935 

7.096 

1 

0.83331 

877.175 

617.53 

325.74 

1028.355 

0.358136 

0.4 

-0,0114 

0.363839 

0.02973227 

620 

Ci.858 

0.32255 

0.2304 

7.044 

1 

0.81172 

877.175 

628.07 

325.74 

1053.211 

0.440925 

0,3 

0.098499 

0.391676 

0,02789959 

630 

Ci.871 

0.38168 

0.2726 

6.985 

1 

0.78668 

877.175 

638.63 

325.74 

1078.656 

0.542928 

0.3 

0.215018 

0.435419 

0,02548809 

640 

0.885 

0.44912 

0.3208 

6.918 

1 

0.7576 

877.175 

649.22 

325.74 

1104.807 

0.66985 

0.3 

0.367087 

0.486306 

0.023222 

The  saturation  temperature  occurs  at  610  K.  This  is  determined  because  at  this 
point  the  difference  between  the  two  values  for  B  is  only  0.01 14.  The  liquid  length  value 
is  read  at  this  point  and  is  29.73  mm  for  hexadecane  at  an  ambient  temperature  of  850  K 
and  amlDient  gas  with  a  density  of  30.2  kg/m^. 

While  performing  these  calculations  the  difference  in  the  B  values  drives  the 
accuracy  of  the  liquid  length.  When  these  values  differ  by  hundredths  the  liquid  length  is 
only  accurate  to  1  mm.  For  this  reason  the  computer  program  should  be  wrote  so  that  the 
B  values  are  within  at  least  thousandths  of  each  other  to  ensure  the  most  accurate  liquid 
lengths. 

This  process  is  repeated  for  each  of  the  five  ambient  temperatures  and  the  five 
ambient  densities  for  the  three  fuels.  In  order  to  check  this  process  hexadecane  is  done 
first;.  The  data  collected  from  the  excel  sheet  is  then  compared  to  the  values  that  are  read 
off  a  chart  in  the  SAE  paper  wrote  by  Siebers.  Table  2  compares  the  liquid  lengths  for 
hexadecane: 


Table  2:  Liquid  length  values  compared  to  previous  research  done  for  hexadecane  at  the 
stated  conditions  _ 


Liquid  Lengths  [mm] 

Ambient 

Temperature 

Density  [kg/m''3] 

3.6 

Seiber's  Calculated 

Seiber's 

7.3 

Calculated 

14.8 

Seiber's 

Calculated 

30.2 

Seiber's  Calculated 

Seiber's 

59 

Calculated 

_ 

Researc;h 

Results 

Research 

Results 

Research 

Results 

Research 

Results 

Research 

Results 

. 

. 

86 

89.89 

63 

66.31 

48 

47.16 

37 

37.01 

850 

88 

90.47 

60 

61.49 

41 

42.07 

29 

29.72 

21 

20.73 

72 

71.43 

47 

48.14 

32 

31,71 

21 

21.38 

15 

15.67 

1150 

61 

62.79 

40 

39.48 

28 

25.88 

17 

17.21 

12 

12.35 

1300 

55 

54.46 

35 

35.40 

23 

23.08 

14 

15,18 

10 

9.76 

]'X  can  be  seen  from  this  table  that  the  values  compare  relatively  well.  Any 
discrepancies  can  be  explained  by  the  difference  in  ambient  conditions  as  well  as  the  fact 
that  Seiber’s  values  for  the  liquid  length  were  estimated  by  reading  them  off  a  chart. 

The  results  for  the  three  fuels  are  summarized  in  the  following  tables  and  figures. 
It  is  seen  that  all  the  fuels  follow  the  same  basic  trends  and  are  all  reasonable. 

Dodecane: 


Table  3 :  Liquid  length  values  for  dodecane  at  the  stated  conditions 


Figure  1 :  The  liquid  length  of  dodecane  at  varying  ambient  temperatures  and  densities 
Tetradeicane: 
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Figure  2:  The  liquid  length  of  tetradecane  at  varying  ambient  temperatures  and  densities 
Hexadecane: 


Table  5 :  Liquid  length  values  for  hexadecane  at  the  stated  conditions 


Liquid  Length  [mm] 

Temperature 

Density  [kg/m''3] 

IKl 

3.6 

7.3 

14.8 

30.2 

59 

700 
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89.89 

66.31 

47.16 

37.01 
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90.47 
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20.73 
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71.43 

48.14 

31.71 
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Figure  3:  The  liquid  length  of  hexadecane  at  varying  ambient  temperatures  and  densities 


Conclusion 

I  t  is  determined  that  this  method  of  determining  the  liquid  length  in  the  three 
diesel  fuel  surrogates  is  accurate.  By  using  the  equations  previously  found  as  well  as 
thermod3mamic  properties  the  liquid  lengths  can  successfully  be  determined.  This  is 
coni'irmed  by  comparing  the  data  with  the  results  previously  obtained  by  Siebers. 

(Tnce  the  process  is  confirmed  the  procedure  is  repeated  for  the  two  remaining 
fuels.  Some  methods  for  finding  the  properties  of  the  fuel  in  the  combustion  chamber  are 
then  simplified  using  best  fit  curves.  These  simplifications  were  then  tested  by  running 
the  calculations  using  the  new  equations  and  comparing  them  to  the  original  data.  It  was 
found  that  the  differences  in  the  liquid  lengths  caused  by  these  simplifications  are 
negligible. 

In  the  end  it  is  determined  that  the  procedures  outlined  in  this  paper  are  suitable 
for  use  in  the  creation  of  a  computer  program.  The  program  should  successfully  be  able 
to  predict  the  liquid  length  in  diesel  fuels. 
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